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Abstract. We consider the Singlet Majorana fermion dark matter model, in which the 
standard model particles interact with the dark sector through the mixing of a singlet scalar 
and the Higgs boson. In this model both the dark matter and the singlet scalar carry lepton 
number, the latter being a bilepton. We suppose the existence of a Z 2 symmetry, remnant of 
some high energy local symmetry breaking, that stabilizes the Majorana fermion. We analyzed 
the parameter space of this model and found that the lepton number symmetry breaking scale, 
drove by the singlet scalar, is constrained to be within hundreds to thousands of GeV, so as to 
conform with the observed dark matter relic density. Finally, we address the direct detection 
and invisible Higgs decay complementarity, confronting our model with recent LUX and LHC 
constraints, as well as XENON1T prospects. 


1. Introduction 

Since almost a century ago, a huge ammount of evidence has been accumulated through 
astrophysics observations 01 El 13 0, 5, 6], indicating that the majority of matter in the 
Universe is of unknown origin. Such observations, however, refer specifically to the 
gravitational response of ordinary (or baryonic) matter to this so called Dark Matter (DM), 
scarcely showing other properties it should comply with, except that it has to be effectively 
electrically neutral [|J non-baryonic, cosmologically stable and cold at recombination. These 
are important clues in order to start building a model to describe the DM. Particle Physics 
is strongly appealing in offering a DM candidate, since many motivating extensions of the 
standard model of electroweak interactions (SM) possess new fields and symmetries, which 
could suitably provide one (or more) stable particle(s) in their spectrum. 

The first observational requirement that a DM model has to fulfill is a a set of interactions 
that can account for the relic density inferred from the power spectrum of the cosmic 
microwave background radiation (CMB). As it says nothing about the particle nature of 
the DM candidate, such as its mass, spin or charges, we have to count on the possibility 
of detecting a signal from a DM particle and gathering more information in order to put 

I Millicharged DM candidates exist in the literature which are perfectly acceptable though mmm 
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this particle under siege. Such a signal should be identified in underground detectors (direct 
detection search), in sattelites or ground-based telescopes (indirect detection search) as well 
as in colliders or accelerators (collider search). There are many experiments searching for 
DM candidates, especially the weakly interacting massive particles (WIMPs). Currently, no 
unambiguous positive DM signal was announced by the experiments, but limits were put on 
its scattering cross-section off nucleon and production rate. 

Any viable particle physics model that addresses the DM issue must comply with these 
bounds and may be constructed in an top-down or a bottom-up way. The top-down models 
intend to explain open particle physics problems in a fundamental way while having a DM 
candidate, such as supersymmetric and extra-dimension models. The bottom-up models 
intend to minimally explain the DM physics, with just a few free parameters and strong 
predictions, and are more easily ruled out by the experiments. This strategy, that we will 
adopt here, allows us to identify what initial hypothesis can be modified in order to agree with 
the experimental limits. 

Models in which the interaction between DM and standard particles is realized through 
scalar exchanges, the Higgs portal, are phenomenologically interesting since the discovery of 
the Higgs-like scalar at the LHC [[TOlCElll. As its couplings and branching ratios are being 
measured, the Higgs physics can constrain the parameter space of such DM models. Here, we 
focus on the fermionic DM in a minimal Higgs portal model. Fermionic DM particles may be 
Majorana or Dirac particles, differing in their relic abundance which is twice larger for Dirac 
fermions, to account for the fact there are two distinct particles annihilating liT2ll . Also, the 
scattering cross-section off nucleon of a Majorana fermion is twice larger relative to a Dirac 
fermion [13]. In many respects there is no significant difference between Dirac and Majorana 
fermion DM candidates though, as pointed out in previous works ITT41 IT2l IT5l 16, 17]. 

Singlet Dirac DM candidates were vastly explored in the context of Higgs portal models 
after the Higgs discovery [ 1^ [19l l20l |2T] |22l [14] . It is known that a lighter than 100 GeV Dirac 
DM is excluded by relic, direct detection and collider constraints, except at the resonance 
region [ 1231 l24ll. Although less considered, a Majorana DM in renormalizable Higgs portal 
models were also pointed out as viable DM candidates 123 (16, 15]. 

The effective field theory (EFT) approach of the Majorana fermion DM was considered 
in [l26l 12,15]. Within EFT approach, the fermionic DM Higgs portal CP-conserving scenario 
is excluded by direct detection search for the 60 GeV - 2 TeV DM mass range. Nevertheless, 
it remains viable if we consider UV completions such as a resonant or an indirect Higgs 
portal Ifl5l . Since we do not know the mass scales of the possible intermediate particles, 
we can not know if an EFT approach is valid in the context of Higgs Portal DM models. 
Indeed, even though renormalizable models must agree with EFT models in the limit of 
heavy mediators, some degrees of freedom concerning the mediator phenomenology are only 
appreciable in concrete models. The limitations of the EFT approach to Higgs Portal DM 
models, concerning the LHC and LUX constraints ll27ll . were considered in refs. Ifl3l 281. 
Also, the implications of the 125GeV scalar for Higgs Portal DM models were studied in a 
model independent way in ref. |[29ll . considering generic scalar, Majorana fermion and vector 
DM candidates. They found that a non-renormalizable Majorana fermion DM is excluded by 
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XENONIOO, if the Higgs-DM coupling strength is larger than 10 -3 . 

In this work, we consider a sterile Majorana neutral fermion as DM candidate in 
a renormalizable CP-conserving Higgs portal model EOt . under the light of the recent 
experimental limits. We suppose the existence of a discrete Z 2 symmetry responsible for 
the Majorana fermion stability, that is expected as a remnant of some spontaneously broken 
gauge symmetry at a high scale not relevant for the low energy physics. The mediator between 
the Majorana fermion and the standard particles is a complex singlet scalar that develops a 
nonzero vacuum expectation value (vev), thus breaking lepton number symmetry. 

This work is the first consideration of the LUX impact on the parameter space of a 
concrete Majorana DM Higgs portal model. We also consider the XENON1T prospect bound 
EH, the invisible Higgs decay constraint and comment how the extension of the scalar sector 
alters the Higgs self-couplings, showing that its future measurement will put important bounds 
on our parameter space. 

This work is organized as follows: in the next section we will present our model. 
In section [3] we compute the relic density and scattering cross section off nucleon for 
the Majorana fermion, in order to constrain our parameter space and discuss how the 
free parameters influence this observables. We also comment the impact of the future 
measurements of the Higgs self-couplings on the free parameters and the possibility of having 
a dark radiation candidate in our scenario. Our main results are presented in section [4| where 
we update the parameter space taking into account the invisible Higgs decay, the LUX and the 
XENON1T bounds, considering separately the scalar extension and the mass spectrum parts. 
Finally, in section [5] we give our conclusions. 

2. The model 

The model we consider here consists of a singlet Majorana right-handed neutral fermion, 
N r , as a DM candidate, whose interaction with the SM particles is through the Higgs portal 
by including a singlet scalar in the model, a, which mix to the SM Higgs doublet, o. This 
Majorana fermion is sterile and stable, a requisite for a DM candidate, when non-trivially 
transforming under a Z 2 symmetry, which forbids it to mix to the light neutrinos and makes it 
interact only with the singlet scalar field. 

We built our Lagrangian assuming lepton number conservation. Considering that N R 
carries one unit of lepton number, no explicit mass term for it is present in the lagrangian. 
Its mass is generated from an Yukawa term involving the X 2 -even singlet bilepton, <7, when it 
acquires a nonzero vev and break the lepton number symmetry. This assumptions reduce the 
number of free parameters in our model, in comparison with the Majorana fermion models 
we referenced. 

The most general renormalizable Lagrangian to be added to the SM, under the above 
assumptions, is given by 

£ D £kin(N R , o) - \ n (N*N r o + 'N r N r o*) - v((f>, a), 


( 1 ) 
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where Cki n (N R , a) is the kinetic term for the new singlet fields and the scalar potential is 

V(f, a) = + A <t>((f) ] (f >) 2 + /i 2 crV + A (T (crV) 2 

+ A^(0 t 0)(crV). (2) 

After spontaneous breaking of electroweak and lepton number symmetries, when (<£>°) = 
v^/s/2 and (o) = v a js/2, the above potential leads to the minimum conditions: 

/'o + V( a2 )o 


and 


v\ = 


- 2 = 


Aw 


Wct + -^0cr(0 2 )o 

A. 


(3) 


(4) 


The physical spectrum is then obtained by shifting the scalar fields to the physical fields, that 
can be written in the unitary gauge as, 



and 


R& T iJ 


a = 


V2 


(5) 


( 6 ) 


Observe that the breaking of lepton number symmetry implies the existence of a Majoron 
in the spectrum, J, which was shown to be a phenomenologically safe Goldstone boson Poll . 
The mass matrix of the massive scalar fields is given by 

/\ < i /: i 

(7) 


j^j-2 _ j 2AA fy G v^v a 

A^Nj^J^cr 2\ rT V„ 


'cr'-'a 


After diagonalization, we have the physical states 

H c a f?0 T SftRfj 
S s a Rfp •' * c a R a , 

where c a and s a are the cosine and sine of the mixing angle a, determined by 


tan 2a = 


\„v 2 ' 


( 8 ) 


(9) 


^4> v <t> ''v^o 

Since the discovery of a scalar resonance in the CMS and ATLAS experiments at 
LHC Mm, with mass Mjj = 125 GeV, and whose interactions indicate it is most probably 
the SM Higgs boson, we can be sure that this scalar cannot be a singlet under SU(2) L . In 
our case, it has to come from the doublet-like scalar, R. 0 , implying that the mixing angle in 
eq. <[8|) has to be small. This is further corroborated by the null search results for a lighter than 
125 GeV scalar, once the coupling of SM fields and the scalar singlet are suppressed by s a . 
We then identify the mass eigenstates by 


M h = 2(A ^c a + A a v a s a + A^tvs^c; 


2 i) 


( 10 ) 
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for the Higgs scalar and 

Mj = 2(A a v 2 a c 2 a + \<t>vls 2 a - A ^v a s 2 a cl). 

for the singlet scalar. 

In terms of the mass eigenvalues, the mixing angle is defined by 

sin 2a = 2X ^ V ^ 


(ID 


m h~ m s' 

and the mixing coupling is given by 


A (j)(j 


M 


H 


M 2 


V4>v a 


( 12 ) 


(13) 


As a consequence of the mixing (eq. [8]), all the standard couplings to the Higgs boson will 
be rescaled in our model by c a . It allows us to probe this model at colliders, looking 
for deviations from the standard interactions with the Higgs boson. For further numerical 
convenience, we can solve the eqs. ( flQ| ) and ( [IT] ) for the quartic couplings, giving 


A fh 


clM 2 + slMl 


A CT = 


clM 2 + s 2 a M* 


(14) 


H 2v ° 

The tree level stability of the scalar potential is guaranteed once A^ > 0, A^. > 0 and 
—2y / A^V < A fa < 2^JX r f ) X r j, which we will use throughout this work[§] 

To complete the spectrum we write down the Majorana fermion mass after the 
spontaneous breaking of lepton number, 

Mn = \/2A n v <j- (15) 

The non-standard interactions introduced by this extension of SM are presented in figure [TJ 



H,S 


H,S,J 


H,S 

/ 

' H,S 
H, S,J 


H , S,J 



H, S,J 


Figure 1. New interactions of the SM extension we are considering. 

§ The stability of a somewhat similar model at one loop level, where the scalar is real and the DM is a Dirac 
fermion, was considered in Ref. [24] while, in Ref. ED, the DM issue was also analyzed along with the 
possibility of a first order phase transition. This last feature not reproducible by our model. 
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M n (GeV) 

M s (GeV) 

v a (GeV) 

a 

[1,20001 

[0.1,20001 

[100,2000] 

[-0.5735,0.5735] 


Table 1 . Ranges for the free parameters used in our numerical scans. 


In summary, from the twelve free parameters, Ajv, /i|, A^, X 0rJ , Mn, Ms, Mh, v^, v a 

and a, when we consider the six constraints, eqs. ([3j|4}|9}[T0}[TT]and[T5]), along with the mea¬ 
sured Higgs mass Mh = 125 GeV and the electroweak scale, v 0 = 246.22 GeV, we may 
choose to keep the remaining independent set of parameters as, M N , Ms, v a and a. In what 
follows, we discuss the role of each of these four free parameters in the observables. 


3. Constraining the parameter space 


Our parameter space is constrained by many requirements. First of all, the Majorana fermion 
must be sufficiently abundant in order to be a viable DM candidate. Here, we use the latest 
results of the Planck telescope l32ll . We computed the relic density and the scattering cross 
section off nucleon (Direct Detection) for the Majorana fermion, using the numerical package 
micrOMEGAs Il33l and compared it with the LUX bounds 1271 . the strongest we have to 
date, as well as with the XENON IT prospects [3_1 j. In this section, we will study how the free 
parameters affect these observables. 

In order to carry out numerical scans on the parameter space, we randomly sampled all 
the free parameters in the ranges shown in table [T] In ref. 041 . we find that the unitarity and 
perturbativity bounds together require Mjy, Ms < 3 TeV and v <7 < 2.4 TeV[U] 

Throughout this work, all the points in the scans will satisfy these conditions: 






a goes from -0.5735 to 0.5735, to be consistent with the LHC bound for the mixing 

(| cos (a) | > 0.84 051 1: 


A^ > 0, X a > 0 (automatic, see eq. 14) and — < X^ < 2y / A^A^ for a potential 

bounded from below; 


• Perturbativity: X 0a , X 0 , X a < 4tt. 


and keep m H = 125 GeV and v a = 246GeV. 

It is instructive to emphasize that this parametrization allows us to consider the two 
possible scalar mass hierarchies: for a positive (negative) A^ CT , a > 0 (a < 0) ensures that 
Ms < Mh and a < 0 (a > 0) ensures that Ms > Mh, as we can see in figure |2} In 
this figure, we see the projection of all the parameter space (with the free parameters varying 
according to the Table [TJ onto the (X 0(J ,a) plane. 

An important feature of this Higgs Portal model is that the ongoing characterization of 
the discovered Higgs boson can provide valuable constraints on the free parameters. One 
example of this is the future measurement of the Higgs sell-couplings If361 :37. 38, 391. which 


|| Although these results are based on the assumption that the Majorana mass terms for fermionic DM are absent 
at tree level, we can naively assume they may be used in our case. 
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Figure 2. Projections of the parameter space onto the plane, for a lighter (left) and 

heavier (right) than Higgs singlet scalar. Note that the former case is possible only for a very 
weak mixing coupling, there is no point for | | much greater than 0.2 under the conditions 

we are considering and the chose range of free parameters. 


we intend to examine next. The 3-Higgs and 4-Higgs self-couplings in our model are given 


by 



(16) 


and 



(17) 


In figure [3} we see the projections of the parameter space onto the (A hhh, ot) (left) and 
(M s , A hhhh) (right) planes. For the 3-H self coupling, we see that deviations from the SM 
value (Xhhh = 31.73 GeV) give us information about the intensity of the scalar mixing. 
Following [35], we see in the 3-Higgs self-coupling projection the bands corresponding to 
the future collider sensitivities: deviations from the SM predicted value of ±50% (TLEP500 
/ HL-LHC-min) are in blue, ±30% (TLEP240 / CEPC HL-LHC-max) in magenta, ±13% 
(ILC) in green and ±5% (hadronlOOTeV) in yellow. For the 4-Higgs self coupling, we see 
that significant deviations from the SM value (.\ HHHH = 0.032) will be able to establish 
the scalar mass hierarchy in our scenario. With this picture in mind we can follow with the 
analysis of the DM constraints on the parameter space. 

3.1. Relic Density 

In this work we are assuming dark matter thermal production, thus the Majorana fermion relic 
density is determined by calculating the thermal averaged annihilation cross section for the 
processes shown in figure |4j 

As the dark matter in our scenario interacts just by scalar exchanges, a characteristic 
behavior of the relic density is one sharply peaked region at Mm ~ Mh /2 ~ 63 GeV, 
corresponding to the Higgs production and depending on the scalar mixture, and another one 
at Mjm ~ Ms/ 2, corresponding to the singlet scalar production. As an asymptotic behavior, 
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Figure 3. 3-Higgs (left) and 4-Higgs (right) self-couplings as functions of the mixing angle 
and singlet scalar mass, respectively. Note that deviations from the standard values give us 
information about the intensity of the scalar mixing as well as the scalar mass hierarchy, 
if Ms is smaller or higher than 125 GeV as we emphasize at the right panel. In the left 
panel, the coulored bands are prospect deviations of ±50% (TLEP500 / HL-LHC-min), in 
blue, ±30% (TLEP240 / CEPC HL-LHC-max), in magenta, ±13% (ILC), in green, and ±5% 
(hadronlOOTeV), in yellow, from the SM predicted value (black line). 



Figure 4. Annihilation channels relevant for the Majorana fermion relic density. 



Mw(GeV) 

Figure 5. Relic density as a function of the Majorana fermion mass. The green interval 
represents the Planck bound on the relic density, for any DM mass. The purple points overdose 
the Universe and will not be considered in our analysis. The gray points indicate regions of 
our parameter space in wich the Majorana fermion is not sufficient to account for the total DM 
content of the Universe. 


we expect a gradual decrease in the relic density, as more channels turn out kinetically allowed 
for more massive Majorana fermions. 

In figure[5j we show the projection of the parameter space onto the (M#, Llh 2 ) plane. Our 
parameter scan shows that there is no set of free parameters allowing for sufficiently abundant 
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Majorana DM lighter than few tens of GeV, in agreement with ll26ll . Considering that the 
lepton number symmetry was spontaneously broken at some scale between 100-2000 GeV, 
the singlet scalar mass from 0.1-2000 GeV and a mixing angle ensuring | cos(a)| > 0.84, we 
can have a Majorana fermion as a DM candidate for masses from 30 GeV to 2 TeV. 

The resonance due to the Higgs exchange, at M N ~ 62.5 GeV, is clear in these plots, 
but if it is going to reach the observed relic density depends mainly on the mixing angle 
value. Stronger the scalar mixing, more intense the annihilation in standard particles and 
consequently smaller the relic density for a Majorana fermion at such mass scale. In figure[6] 
we show two particular values of Ms that set different singlet scalar resonances. 


iV=2TeV, M s =100GeV iv=2TeV, M s =1000GeV 




Figure 6. Singlet scalar resonances, for all values of the mixing angle. Notice that the Planck- 
allowed points (green) are always near the scalar resonances and where the abundance curve 
falls down. It is a characteristic behavior of this model, as we will see in more detail in what 
follows. Same color code of figure[5] 

Let us now see how the lepton number symmetry breaking scale, v a , changes the relic 
density. In figure [7} we see that there is a global increase in the relic density according to the 
increase of v a , due mainly to the interaction of the Majorana fermion with H, S and J, which 
is inversely proportional to v a . In view of this, we have an upper and a lower allowed scale for 
the global symmetry breaking to occur in this DM scenario. The plots in figure [7] comprehend 
all the Ms values, then all the singlet scalar resonances. As we are taking singlet scalars up 
to 2 TeV, the resonances due to their productions go up to 1 TeV. 

We can now look at the parameter space under the light of direct detection experiments, 
which we do next. 


3.2. Direct Detection 


The interaction between the Majorana fermion and the quarks are realized via t-channel scalar 
exchanges. So, the scattering cross section is always spin-independent (SI) and is given by 

4M^ 1/2 / /l^2\-l/2 


&NNqq ~ 


\M-NNqqY 


1287 t 2 EL 


1 - 


EL 


-> 4 Ml 

EL 


(18) 


The scattering amplitude is the sum of each scalar contribution: 


M 


NNqq 


= M h + Ms ex i 


MqMN 

V<j>V (7 


SfyCr 


1 


Ml 


Ml 


(19) 
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iV=100GeV 
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M N (GeV) 


1 5 10 50 100 
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Figure 7. Influence of the lepton number symmetry breaking scale on the relic density. Same 
color code of figure [5] 


where we have taken the low energy limit in which the masses of the virtual particles are 
much greater than their momenta. This minus sign comes from the diagonalization of the 
scalar mass eigenstates (eq.[8]). 

We present a scatter plot for the scattering cross section off nucleon in figure [8j that 
shows the well known interference effect in the amplitude near the Higgs mass value, as can 
be expected from eq. [19} This is why the case of Ms ~ 125 GeV evades the LUX constraint, 
as well as the future XENON IT constraint, for all the free parameter values considered. We 
also see the enhancement of the scattering cross section provided by low scalar masses. 

In figure [9] we show some plots that clarify the effects of the singlet scalar mass on the 
WIMP-nucleon scattering cross section, taking into account the points that are excluded by 
the LUX results (red), and the points within the Planck interval (gree) for the relic density, as 
well as the XENON1T exclusion prospect (dark red). 

Although there exist the possibility of a lighter than Higgs singlet scalar in this 
model IflOl . this case is quite ruled out by LUX in our DM scenario, as we can see in figure [9] 
(left). Henceforth, we will focus on the case of a singlet scalar heavier than the Higgs boson. 
As we have already remarked, the case in which the scalars are quite degenerate leads to 
a cancellation effect in the scattering amplitude, evading the LUX constraint and even the 
XENON1T future constraint (figure [9} middle). Notice that XENON1T will be restrictive 
even for a TeV scale singlet scalar (figure [9} right). 

The inviability of very light singlet scalars in this model has an interesting consequence. 
As we have seen, there exists in our spectrum a Majoron, J, that is a massless pseudoscalar. 
Cosmologically, this particle will naturally contribute to the radiation energy density as long 
as its interactions make it decouple from the thermal bath. As there are evidences for an 
unknown contribution for the relativistic degrees of freedom at recombination, the so called 
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I 

i 

S3 


1 5 10 50 100 500 1000 

M s ( GeV) 

Figure 8. Scattering cross section off nucleon as a function of the singlet scalar mass. The 
red points are already excluded by the LUX results. The dark red points will be excluded by 
XENON IT if no signal is confirmed in the next few years. Notice that the case of degenerate 
scalars can evades the direct detection bounds for all the free parameter values and very light 
singlet scalars are not favored since they enhance the scattering cross section. 



dark radiation, we have checked the possibility of J address this question in our dark matter 
scenario. 

The Majoron can be a radiation at the recombination epoch if it decouples just before the 
muon annihilation [4JJ and in that case, roughly speaking, we must have 

™ ~ 10 GeV. (20) 

v \^<f)(T I 

As we are taking only perturbative mixing coupling values, A ^ < 47t and so M s < 
35 GeV. As we can infer from figure [8j this possibility is quite excluded by direct detection 
bounds. 


iV=2TeV, M s =lOGeV v a =2TeV, M s =lOOGeV v„= 2TeV, M s =lTeV 



M n ( GeV) M n ( GeV) M N ( GeV) 


Figure 9. Scattering cross section as a function of the Majorana fermion for all values of the 
mixing angle and certain values of the singlet scalar mass. We see that a lighter than Higgs 
singlet scalar is quite ruled out by LUX (left). The cases of a quite degenerate (middle) and 
a TeV scale (right) singlet scalar can evade the direct detection bounds, but will be strongly 
constrained by the XENON1T results. 


Finally, in figure [10] we see that we can evade the direct detection bounds for a small 
enough mixing angle (left), as the interacion of the Majorana fermion with the standard 
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particles is drove by a Higgs portal. This was also observed in ref. ll2Tll . Again, in the 
scalar degenerate case, when Ms ~ 125 GeV, we can also evade the direct detection bounds, 
but the XENON IT future results should constrain strongly this case as well. In these plots, 
such a scenario has Planck-allowed points at M N ~ 1.2 TeV. All this analysis is a fair 
update concerning the issue of DM relic density and direct detection. Once we have a lot 
of information since the Higgs boson discovery we can use it to further improve our results on 
the parameter space siege for this model, the so called collider complementarity. We perform 
this task in the following section. 


tV=2TeV, te-O.Ol tv=2TeV, a=0.2 i/„=2TeV, a=0.5 



10 20 50 100 200 500 10002000 10 20 50 100 200 500 10002000 10 20 50 100 200 500 10002000 

M n ( GeV) M w ( GeV) M N ( GeV) 

Figure 10. Scattering cross section as a function of the Majorana fermion for certain values of 
mixing angle. The Planck-allowed points are for a heavier than Higgs singlet scalar. 

In summary, we have restricted our parameter space by considering only coupling 
constants within the perturbative regime and ensuring a potential bounded from below. 

We have considered mixing angles compatible with the maximum mixing allowed for the 
discovered Higgs boson and examined the impact of the future measurements of the Higgs 
self-couplings on our model. The values for v a were limited from thousands of GeV to few 
TeV in order to furnish the observed relic density and we have seen that our model is not 
capable of explaining the existence of a very light singlet scalar if we consider it related to the 
dark matter physics, through the direct detection bounds. Now, we are ready to look for viable 
regions of our parameter space taking the experimental bounds in a complementary way. 


4. Complementarity bounds 


This model predict new decay channels for the Higgs boson, that are so far invisible for the 
LHC. The new decay rates are given by: 


r(H JJ) = 

T(H -)• NN) = -L- fljgjgg ( 






T(H SS) = 


Mi ££ 

327T 


1 + 


2 M% 
Mi 


c a + s a ^ 


1 4 Mg 

!- m# 


Considering the Higgs decay rate in SM as T SM ~ 4 MeV [142], the upper bound on the 
branching ratio for the invisible Higgs decay (IHD) ll43Tl is 


TDinV _ _ / A A 

H - clT SM + p inv ' 5 


( 21 ) 
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< 2.67 x 10 ~ 3 GeV. 


( 22 ) 


As we focus in the heavier than Higgs singlet scalar case, there are only two possible 
hierarchies: M H < 2 M s , 2 M N and 2 M N < M H < 2 M s . For T inv = Y(H —y JJ), we have 

t9a ' 1 A in “ 7 (23) 


and for F' 


— < 1.4 x 10“ 7 
v 2 

= T(H -> JJ) + T(H -> NN ), 
Ml A 4Ml \ 3/2 


* 9 , 


1 + 2 


l N 

Ml 


1 


Ml) 


< 1.4 x 10 


-7 


(24) 


As the Z boson do not decay at tree level in any singlet scalar in this model, the constraints 
from invisible Z decay is not relevant here. 

In order to explore the regions of our parameter space allowed by the Planck, LUX, 
XENON1T and IHD constraints, we will choose values of the free parameters based on the 
study of the previous section, instead of arbitrary ones. We do this next by fixing first pairs of 
free parameters that parametrize our scalar sector extension, v a and a, and then pairs of mass 
values of the new particles of our spectrum, Mjv and Ms- 

For all the following free parameter choices, we have checked that the coupling constants 
Atv is always smaller than unity. Also, all the viable parameter space regions are consistent 
with X 0a < 1 because of the direct detection bounds. 

In what follows, we will identify the regions excluded by the invisible Higgs decay 
bounds. The regions filled with blue straight lines are in the hierarchy M H < 2M n . 2 Ms 
and do not obey the eq. [23} The regions filled with orange straight lines are in the hierarchy 
2 Mjv < Mh < 2Ms and do not obey the eq. [24} Consistently with our previous figures, the 
green lines are composed of points within the Planck interval giving a sufficiently abundant 
Majorana fermion, while the gray regions are composed by the points that give underabundant 
Majorana fermion, indicating the regions of the parameter space in which the Majorana 
fermion could be only a fraction of the total DM content of the Universe. Although we do not 
consider another DM candidate or non-thermal production that validate these gray points, we 
keep it for completeness. Regarding the direct detection search, the red regions are excluded 
by LUX and the dashed dark red lines will extend this exclusion if XENON IT do not confirm 
any DM signal in a close future. 


4.1. v a and a fixed 

In this subsection, we fix a set of values for the new vev, v a , and the mixing angle, a, that 
parametrize the scalar extension of this model. Based on our previous results, lustrated by 
figure [7} we will study how the lepton number breaking scale changes the allowed regions 
by choosing a vev v a = 500 GeV, 1 TeV and 2 TeV. The role of the mixing angle will be 
considered by taking it as a = 0.01, 0.2 and 0.5. 
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Figure 11. Free parameter slices (Mn,Ms) for fixed v a values (upper panel) and fixed a 
values (lower panel). The green region are in agreement with Planck while the gray one give 
underabundant Majorana fermion. The red region are already excluded by LUX and the dark 
red dashed lines shows how XENON1T should expand this exclusion if do not observe a DM 
signal. The blue (orange) straight lines show the region excluded by IHD into two majorons 
(two majorons and two Majorana fermions). 


In figure [TTJ we display the free parameter slice {M N , M s ) for different pairs of a and 
v a . All the Planck-allowed points (green) are near the scalar resonance, in agreement with 
previous works. For strong enough scalar mixing (a = 0.2 and a = 0.5 in this figure) we see 
Planck-allowed points for M N ~ M H /2 = 62.5 GeV, corresponding to the Higgs resonance. 
As we are varying M s up to 2 TeV, in the three first plots (for v a = 2 TeV) the Planck- 
allowed points for M N < 1 TeV correspond to the singlet scalar resonances and the ones for 
Mm > 1 TeV, to the falls of the abundance curve (see figure |6j). 
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For a = 0.01 and a = 0.2, the direct detection search do not exclude any Planck- 
allowed point, but for a = 0.5 it offers a strong limit. We see that > 500 GeV is 
excluded by LUX for a = 0.5 if M s > 500 GeV and if XENON1T do not confirm any signal, 
this case becomes practically excluded. For a TeV scale vev, all the points are allowed by 
the IHD, but v a = 500 GeV (left lower panel) is completely excluded by IHD constraint, 
considering a conservative mixing angle of 0.2. For a v a = 1 TeV (right lower panel), we 
see that Mjy > 500 GeV is also excluded by LUX for a = 0.5 if Ms > 500 GeV but in 
comparison with the case “v a = 1 TeV, a = 0.5” (right upper plot), we have more Planck- 
allowed points in the LUX-allowed region. 

In the limit of null mixture (a —> 0) and very high energy breaking scale of the symmetry 
we are supposing (v a —> oo), the scalar sector of our model becomes standard, with just one 
scalar coming from a doublet, no Majoron in the spectrum and without the possibility of 
provides mass for our dark matter candidate. As we have seen in figure |TT| our parameter 
space is favored by small deviations of the scalar sector, it means, for smaller scalar mixing 
and very high energy physics (left upper panel), in view of the current collider and direct 
detection searches. We have seen that the model we consider here is very predictive and in the 
next few decades, all this scenario will be tested by the experiments. It is constructive to keep 
in mind, however, that the confirmations of a possible DM particle signal and a possible new 
scalar discovery at colliders or accelerators are very connected things for Higgs portal models 
and their immersion in more complex contexts may offer important hints. 

4.2. Mn and Ms fixed 

Now we study the allowed parameter space for fixed values of M N and M s . Here we consider 
the lepton number breaking scale, v a , in a larger window, from 100 GeV to 10 TeV. The TeV 
scale is a natural scale in the search for new physics, so we study separately the case of 
Mn = 1 TeV and Ms = 1 TeV. For a M^ = 1 TeV, we show how the allowed region depends 
on M s by fixing it at 130 GeV (quite degenerate case), 500 GeV (intermadiate case) and 1 TeV 
(TeV scale physics). On the other hand, for Ms = 1 TeV, we consider Majorana fermions 
with masses of 60 GeV (Higgs resonance case) and 500 GeV (singlet scalar resonance case). 

In figure |T2j we show the free parameter slices (a, v a ). In the upper and middle panels, 
we have the case of Mjv = 1 TeV. The Majorana fermion will be sufficiently abundant only 
if the lepton number symmetry breaking is near 2 TeV, regardless the singlet scalar mass, 
implying that its interaction with the scalars would be weak (see eq. [15] ). Lower breaking 
scales would need another DM candidate or a non-thermal DM production, to be constrained 
by IHD, LUX and XENON IT. As we already pointed out, the quite degenerate case (left 
upper panel) can evade direct detection and as we may expect, the LHC provides the strongest 
current bound, but cannot constrain the scalar mixing of the Higgs with a scalar that develops 
vev too early, at TeV scale. The XENON1T will allow scalar mixings up to 0.4 in that case. 
For heavier singlet scalars, the current limits for the mixing are a < 0.3 for M s = 500 GeV 
(right upper panel) and a < 0.2 for Ms = 1 TeV (middle panel). The XENON1T future limits 
will restrict much more the mixture, up to a < 0.05, in which case non-standard interactions 


Majorana Dark Matter in Minimal Higgs Portal Models after LUX 


16 


M n = ITeV, M s = 130GeV 


M n = ITeV, M s = 500GeV 

XENON IT \ / 




Figure 12. Free parameter slices (a, v a ) for fixed vev values (upper panel) and fixed mixing 
angle values (lower panel). Same color code of figure [TT] 


will not be detected easily. 

In the lower panel, we can appreciate how the scalar resonance regions of the spectrum 
evade the current experimental bounds. At left, we show the Higgs resonance case (Mn = 
60 GeV), where the LHC bound is stronger than LUX bound but cannot constraint a 
sufficiently abundant Majorana fermion. In this case, the breaking scale should be from 
200 GeV (for small mixings) to 5 TeV (for maximal mixing). Nevertheless, if no signal is 
confirmed by XENON1T, this case will be completely excluded in the next few years. At 
right, we show the singlet scalar resonance (M N = 500 GeV), also not constrained by the 
current bounds but accessible to XENON1T, when the mixing will be restricted to a < 0.1. 
This possibility requires the breaking scale near 5 TeV for a sufficiently abundant Majorana 
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fermion. 

We have seen that even the safer resonance regions will be strongly constrained in the 
near future, and even considering a TeY scale new physics. 


5. Conclusions 


In this work, we updated the Majorana fermion DM in a CP-conserving Higgs portal scenario, 
in view of the LUX, XENON1T and the invisible Higgs decay (IHD) recent results. 

Our dark matter candidate is a sterile neutral Majorana fermion and the Higgs portal 
is opened by a complex scalar field, singlet and neutral under the standard symmetries and 
charges. This singlet scalar develops spontaneously a vacuum expectation value (vev) that 
gives a Majorana mass for the fermion and leaves it stable due to the discrete Z 2 symmetry we 
are supposing, making it a DM candidate. This model is described by only four parameters: 
the Majorana fermion and singlet scalar masses, M N and Ms, the mixing angle between the 
scalars, a, and the lepton number breaking scale, v a . 

Our parameter space is within the sensibility region of the future collider searches in what 
concerns to the 3-Higgs self-coupling constant. It means that the measurement of the Higgs 
self-coupling can shade light on our Majorana fermion DM scenario. Also, strong deviations 
from the SM value of the 4-Higgs self-coupling will give us information about the scalar mass 
hierarchy. All the analysis we have made is within the perturbative regime and is compatible 
with a bounded from below scalar potential. 

We found that v a must be between few hundreds of GeV and few TeV in order to the 
Majorana fermion account for the total amount of DM in the Universe, according to the 
Planck results. Considering that this symmetry breaking took place at some scale between 
100-2000GeV, the singlet scalar mass from 0.1-2000GeV and a mixing angle ensuring 
| cos (a) | > 0.84, we can have a Majorana fermion as a DM candidate for masses from 30GeV 
to 2TeV. The case of a lighter than Higgs singlet scalar in this DM scenario is quite ruled out 
by LUX, making inviable to consider the Majoron as dark radiation. 

With the aim of “split the degeneracy” of this allowed region, we studied the role of 
each parameter to find the very predictive curves allowed by Planck and contrasted by the 
experimental limits. We did it systematically by choosing benchmark values for the free 
parameters in two steps. Firstly we examined the scalar extension parametrization, by fixing 
a at 0.01, 0.2 and 0.5 and v a at 500GeV, ITeV and 2TeV. Secondly, we examined the mass 
spectrum, by fixing M N at 60GeV, 500GeV and ITeV and M s at 130GeV, 500GeV and ITeV. 

As expected, the allowed regions are near the scalar resonances. We have shown that it is 
possible to evade the LUX and even the XENON1T constraint for a Majorana DM candidate 
in two situations: scalar degenerate case (fig. [9]) and small mixing angles (fig. my 

The case of a relatively low breaking scale (v a = 500GeV) was found totally excluded 


by IHD even for a mixing angle of only 0.2 and TeV scale v a values were found safe (fig. [IT]). 
We have shown that Mm > 500 GeV is excluded by LUX if Ms > 500 GeV in essentially 
two situations: “ v a = 2 TeV,a = 0.5” and “v a = 1 TeV,a = 0.2” (right upper and lower of 


% El respectively). 
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For M n = 1 TeV, v a must be at approximately 2 TeV in order to agree with the Planck 
results and the mixing was found restricted to a < 0.3 (for Ms = 500 GeV) and a < 0.2 
(for M s = 1 TeV) by the LUX results and to ct < 0.05 by the XENON IT future bound 
(upper and middle panels of fig. |T2j). It means that if no signal is confirmed by XENON1T, it 
becomes very difficult to detect some non-standard interaction at colliders, since they should 
be proportional to sin ( a ) 2 . The resonance regions were found not constrained by the current 
experiments, but XENON1T will exclude completely M N ~ 60 GeV and restrict the mixing 
to a < 0.1 for M s = 1 TeV (lower panel of fig.[l2|). 

We have seen how the discoveries of a DM and a new scalar particles are strongly related 
in our Higgs portal scenario. While this minimal scheme remains allowed by the experimental 
bounds, it can be embedded in more fundamental extensions of the SM. In supersymmetric 
theories, stable neutralinos are Majorana fermions that also appear as natural dark matter can¬ 
didates and in some gauge extensions and neutral scalars may be the remnants of a more 
fundamental gauge symmetry breaking, as studied in P4ll . 


Acknowledgments: 

We are very grateful to Clarissa Siqueira, Antonio Oliveira and Farinaldo Queiroz for useful 
discussions. This work was supported by Coordenagao de Aperfeigoamento de Pessoal de 
Nivel Superior - CAPES (MD) and Conselho Nacional de Desenvolvimento Cientifico e Tec- 
nologico - CNPq (CASP,PSRS). 


References 


[1] F. Zwicky, Republication of: The redshift of extragalactic nebulae, General Relativity and Gravitation 41 

(Jan., 2009) 207-224. 

[2] V. C. Rubin and W. K. Ford, Jr., Rotation of the Andromeda Nebula from a Spectroscopic Survey of 

Emission Regions, As trophy s.J. 159 (Feb., 1970) 379. 

[3] M. A. Monroy-Rodriguez and C. Allen, The end of the macho era, revisited: New limits on macho masses 

from halo wide binaries, Astrophys.J. 790 (2014), no. 2 159. 

[4] D. Clowe, M. Bradac, A. H. Gonzalez, M. Markevitch, S. W. Randall, et al., A direct empirical proof of the 

existence of dark matter, Astrophys.J. 648 (2006) L109-L113, |astro-ph/0 6084 07]. 

[5] A. Klypin and F. Prada, Testing gravity with motion of satellites around galaxies: Newtonian gravity 

against Modified Newtonian Dynamics, Astrophys.J. 690 (2009) 1488-1496, [ arXiv : 070 6.3554]. 

[6] S. Dodelson, The Real Problem with MOND, International Journal of Modern Physics D 20 (2011) 

2749-2753, [arXiv: 1112.1320]. 

[7] H. Goldberg and L. J. Hall, A new candidate for dark matter, Physics Letters B 174 (1986), no. 2 151 - 155. 

[8] A. de Rujula, S. L. Glashow, and U. Sarid, Charged dark matter, Nuclear Physics B 333 (Mar., 1990) 

173-194. 

[9] L. Chuzhoy and E. W. Kolb, Reopening the window on charged dark matter , JCAP 0907 (2009) 014, 

[arXiv : 0809.043 6]. 

[10] CMS Collaboration, S. Chatrchyan et al., Observation of a new boson at a mass of 125 GeV with the CMS 

experiment at the LHC, Phys.Lett. B716 (2012) 30-61, [arXiv: 1207.7235]. 

[11] ATLAS Collaboration, G. Aad et al., Observation of a new particle in the search for the Standard Model 

Higgs boson with the ATLAS detector at the LHC, Phys.Lett. B716 (2012) 1-29, [arXiv: 1207. 721 4]. 

[12] M. A. Fedderke, J.-Y. Chen, E. W. Kolb, and L.-T. Wang, The Fermionic Dark Matter Higgs Portal: an 

effective field theory approach, JHEP 1408 (2014) 122, [arXiv: 14 04.22 83]. 



Majorana Dark Matter in Minimal Higgs Portal Models after LUX 


19 


[13] S. Baek, P. Ko, and W.-I. Park, Invisible Higgs Decay Width vs. Dark Matter Direct Detection Cross 

Section in Higgs Portal Dark Matter Models, Phys.Rev. D90 (2014), no. 5 055014, 

[arXiv: 1405.3530]. 

[14] S. Baek, P. Ko, and W.-I. Park, Search for the Higgs portal to a singlet fermionic dark matter at the LHC , 

JHEP 1202 (2012) 047, [arXiv -1112.1847 ]. 

[15] L. Lopez-Honorez, T. Schwetz, and J. Zupan, Higgs portal, fermionic dark matter, and a Standard Model 

like Higgs at 125 GeV , Phys.Lett. B716 (2012) 179-185, [arXiv: 1203.2064]. 

[16] N. Okada and T. Yamada, Simple fermionic dark matter models and Higgs boson couplings , JHEP 1310 

(2013) 017, [arXiv: 1304.2962]. 

[17] T. Boeckel and J. Schaffner-Bielich, Cosmology of fermionic dark matter , Phys.Rev. D76 (2007) 103509, 

[arXiv : 0707.32 60]. 

[18] K. Ghorbani, Fermionic dark matter with pseudo-scalar Yukawa interaction , JCAP 1501 (2015) 015, 

[arXiv : 1408.492 9]. 

[19] Z. Bagherian, M. M. Ettefaghi, Z. Haghgouyan, and R. Moazzemi, A new parameter space study of the 

fermionic cold dark matter model , JCAP 1410 (2014), no. 10 033, [arXiv: 14 0 6.2 92 7]. 

[20] T. Franarin, C. Vasconcellos, and D. Hadjimichef, On the possibility of a 130 GeV gamma-ray line from 

annihilating singlet fermionic dark matter , Astron.Nachr. 06 (2014) 07, [arXiv: 14 04.0 4 0 6]. 

[21] M. Fairbairn and R. Hogan, Singlet Fermionic Dark Matter and the Electroweak Phase Transition , JHEP 

1309 (2013) 022, [arXiv: 1305.3452]. 

[22] C.-K. Chua and R.-C. Hsieh, Study of Dirac fermionic dark matter , Phys.Rev. D88 (2013), no. 3 036011, 

[arXiv: 1305 770 08]. 

[23] Y. G. Kim, K. Y. Fee, and S. Shin, Singlet fermionic dark matter , JHEP 0805 (2008) 100, 

[arXiv :0803.2932[. 

[24] S. Baek, P. Ko, W.-I. Park, and E. Senaha, Vacuum structure and stability of a singlet fermion dark matter 

model with a singlet scalar messenger , JHEP 1211 (2012) 116, [arXiv: 120 9.4163]. 

[25] S. Esch, M. Klasen, and C. E. Yaguna, Detection prospects of singlet fermionic dark matter , Phys.Rev. 

D88 (2013) 075017, [arXiv: 130 8.0 951]. 

[26] S. Matsumoto, S. Mukhopadhyay, and Y.-L. S. Tsai, Singlet Majorana fermion dark matter: a 

comprehensive analysis in effective field theory , JHEP 1410 (2014) 155, [arXiv: 14 07.185 9]. 

[27] LUX Collaboration, D. Akerib et al., First results from the LUX dark matter experiment at the Sanford 

Underground Research Facility , Phys.Rev.Lett. 112 (2014) 091303, [ arXiv: 1310.8214]. 

[28] G. Busoni, A. De Simone, E. Morgante, and A. Riotto, On the Validity of the Effective Field Theory for 

Dark Matter Searches at the LHC , Phys.Lett. B728 (2014) 412-421, [arXiv: 1307.2253]. 

[29] A. Djouadi, O. Lebedev, Y. Mambrini, and J. Quevillon, Implications of LHC searches for Higgs-portal 

dark matter , Phys.Lett. B709 (2012) 65-69, [ arXiv: 1112.32 9 9]. 

[30] C. de S. Pires, F. Queiroz, and P. Rodrigues da Silva, Singlet Majorana fermion dark matter, DAMA, 

CoGeNT, and CDMS-If Phys.Rev. D82 (2010) 105014, [arXiv: 1002.4601]. 

[31] E. Aprile and XENON1T collaboration, The XENON IT Dark Matter Search Experiment , ArXiv e-prints 

(June, 2012) [arXiv: 1206.6288]. 

[32] Planck Collaboration, P. Ade et al., Planck 2013 results. XVI. Cosmological parameters , 

Astron.Astrophys. 571 (2014) A16, [arXiv: 1303.507 6]. 

[33] G. Belanger, F. Boudjema, A. Pukhov, and A. Semenov, micrOMEGAs-3: a program for calculating dark 

matter oservables, Comput.Phys.Commun. 185 (2014) 960-985, [arXiv: 1305.0237]. 

[34] D. G. E. Walker, Unitarity Constraints on Higgs Portals , arXiv:1310.1083 

[35] S. Profumo, M. J. Ramsey-Musolf, C. L. Wainwright, and P. Winslow, Singlet-catalyzed electroweak 

phase transitions and precision Higgs boson studies , Phys.Rev. D91 (2015), no. 3 035018, 

[arXiv :1407.5342[. 

[36] A. J. Barr, M. J. Dolan, C. Englert, D. E. Ferreira de Lima, and M. Spannowsky, Higgs Self-Coupling 

Measurements at a 100 TeV Hadron Collider , JHEP 1502 (2015) 016, [arXiv: 1412.7154]. 

[37] A. Conway, H. Wenzel, R. Lipton, and E. Eichten, Measuring the Higgs Self-Coupling Constant at a 

Multi-TeVMuon Collider , arXiv: 1405.5910, 


Majorana Dark Matter in Minimal Higgs Portal Models after LUX 


20 


[38] ILD Collaboration, J. Tian and K. Fujii, Measurement of Higgs couplings and self-coupling at the ILC , 

PoS EPS-HEP2013 (2013) 316, [arXiv: 1311.652 8]. 

[39] F. Goertz, A. Papaefstathiou, L. L. Yang, and J. Zurita, Measuring the Higgs boson self-coupling at the 

LHC using ratios of cross sections , arXiv:1309.3805, 

[40] J. D. Clarke, R. Foot, and R. R. Volkas, Phenomenology of a very light scalar (100 MeV < rrih < 10 GeV) 

mixing with the SMHiggs , JHEP 1402 (2014) 123, [arXiv: 1310.80 42]. 

[41] S. Weinberg, Goldstone Bosons as Fractional Cosmic Neutrinos , Phys.Rev.Lett. 110 (2013), no. 24 

241301, [arXiv: 1305.1971]. 

[42] CMS Collaboration, V. Khachatryan et al., Precise determination of the mass of the Higgs boson and tests 

of compatibility of its couplings with the standard model predictions using proton collisions at 7 and 8 
TeV, arXiv : 1412.8 662, 

[43] N. Zhou, Z. Khechadoorian, D. Whiteson, and T. M. Tait, Erratum: Bounds on Invisible Higgs boson 

Decays from ttH Production , Phys.Rev.Lett. 113 (2014) 151801, [arXiv: 1408.0011]. 

[44] P. Dong, D. Huong, F. S. Queiroz, and N. Thuy, Phenomenology of the 3-3-1-1 model , Phys.Rev. D90 

(2014), no. 7 075021, [arXiv: 14 05.25 91]. 


